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Abstract. 

We consider the N-soliton solutions in the sine-Gordon model as a N-body problem. This 
leads to a relativistic generalization of the Calogero model first introduced by Ruijsenaars. 
We show that the fundamental Poisson bracket of the Lax matrix is quadratic, and the r- 
matrix is a dynamical one. This is in contrast to the Calogero model where the fundamental 
Poisson bracket of the Lax matrix is linear. 
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1 Introduction. 



We reconsider the relation between the sine-Gordon solitons and the relativistically invariant 
N-body problem introduced by Ruijsenaars. One of the motivations consists in extracting 
coordinates suitable for quantization. Since in the quantum theory, the asymptotic solitons 
are of primary importance [p], § , soliton related coordinates seem appropriate. As illustrated 
by eq.([IT]), the sine-Gordon field is expressed in a simple way in terms of these coordinates. 
This kind of formula could be useful for computing N-soliton contribution to e.g. Casimir 
energy, form factors, etc... 

The Ruijsenaars models H are relativistic generalizations of the Calogero models Q. 
The later are integrable and even exactly solved: the spectrum and more remarkably the 
wave functions are known. Recently their rich algebraic structure began to appear ||. It is 
tempting to believe that a similar structure will still exist in these relativistic models. 

The feature of the Calogero models that we want to extend to the relativistic case is the 
existence of a r-matrix ||. We find that the Poisson bracket for the Lax matrix is quadratic, 
in constrast to the non-relativistic case where it is linear. The existence of this quadratic 
Poisson bracket, eq.(^), is probably a manifestation of the quantum afline symmetry of the 
sine-Gordon field theory 0. Despite the fact that our motivation is the quantum theory, in 
this Letter we restrict ourselves to the classical theory. We illustrate the use of the r-matrix 
by providing simple proofs for the commuting property of the Hamiltonians and for the 
symplectic property for a remarkable coordinate transformation. It is interesting to notice 
that after quantization, the Hamiltonians eq.(2IJ22) will lead to finite difference problems 
instead of differential equations. 

We would like to point out that elliptic generalizations of these models exist. It was 
recently shown independently by Sklyanin, that a quadratic Poisson bracket also exists for 
their Lax matrix IHI. 



2 Tau functions and solitons. 

We first introduce few notations for the sine-Gordon equation and its solutions. Let z± = x±t 
be the light cone coordinates and d± = \{d x ± d t ). The sine-Gordon equation is: 

<9+<9_0 = 2sin(20) (1) 

It is convenient to introduce the two tau functions r±, which satisfy the Hirota equations: 

r±{d-d + T ± ) - (9 + r ± )(5_r±) = tI-t% (2) 

The sine-Gordon field <fi is related to the tau functions by — = exp(— i<p). 

The tau functions of the iV-soliton solutions of the sine-Gordon equation are given by 

asm 

Ti N \z +J z.) =det(l±V) 
with V a. N x N matrix with elements: 



Vij = 2^^-yfx^Cj with X i = a i ex9(2( l i i z+ + nT 1 z-)) (3) 
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The quantities dj and m are the parameters of the solitons: yUj are the rapidities and are 
related to the positions. For the sine-Gordon equation, they satisfy specific reality conditions. 
For solitons or antisolitons, the rapidity fi is real and a is purely imaginary, i.e. a = iee 1 
with e = +1 for a soliton and e = — 1 for an antisoliton. The "breathers" correspond to pairs 
of complex conjugated rapidities (fi,~p) and positions (a, —a). Notice that these conditions 
are preserved by the dynamics. 

The sine-Gordon equation is a Hamiltonian system. The symplectic form is the canonical 
one: 

/+oo 
dx 5ti(x) A S(j)(x) 
-oo 

with tt(x) the momentum conjugated to the field <ft(x). Above, 5 denote the differential 
on the phase space. One can take the restriction of this symplectic form to the iV-soliton 
subspace of the phase space. In the coordinates aj and /ij, we find that this restriction is: 

ii »■ in t^j \/''r • • t<> t'j 

This two-form is non degenerate. It therefore defines a symplectic structure on the restricted 
phase space. The corresponding Poisson brackets are found by inverting the symplectic form 

{a^/ij} = a^jdij (5) 




{a i} aj} = - 

How can Eq.(^j) be proved? First, we consider the one and two soliton solutions. In these 
cases, the computation can be done directly using the formula for the field <p(x) and the 
momentum n(x) = d t (j)(x). We find: 



n 



SG 



restricted 




u (6) 



Thanks to the behaviour of the r-functions at infinity, the prefactor is an irrelevant numerical 
constant. Then, we consider the general case with an arbitrary number of solitons with 
parameters (/ij, a^). Since the dynamics is Hamiltonian, the symplectic form can be computed 
at any time; in particular at t — ► ±oo. In this in- or out-limit, the sine-Gordon field becomes 
asymptotically equal to the sum of the one-soliton solutions with parameters (/z- n , af 1 ) and 



/4 n = /i° ut = ^ 



ar = «i n 

I < I 



Ia«jI> Iwl 

Since asymptotically the solitons decouple, in the symplectic form the crossed terms vanish 
(the overlap integrals are zero). Therefore the symplectic form reduces to the sum of the 
one-soliton expressions, but with the shifted in and out parameters: 

w = y 6 < A 5 _B - y 6a ° ut A gAg 

i=l ®i f^i i=l f^i 
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This is equal to eq.(|]). If breathers are present the proof is identical since we checked eq.(|J) 
up to two solitons. 

As a by product, we see that the transformation from the in-variables (/z™, a™) to the 
out-variables (/x° ut , a° ut ) is symplectic. The classical S-matrix is the generating function for 
this transformation. 

Now that we have specified the Poisson brackets between the parameters (fii, <2j), we can 
compute the Poisson brackets between the matrix elements V^. Remarkably, the latter can 
be written with the help of a <i-matrix depending on the dynamical variables. 

Proposition. Denote by Eij = the canonical basis of the N x N matrices. Put 

V = Eij VijEij. Then: 

{Vi f ^2} =[di2,K]-[d2i,^2] (7) 
where d 12 = <>,, 11 1'-,;, ® E u with, 

dij. k i = ~ [^Z^f] (VikSu + VjtSik + V ik 5 3l + V a 5 jk ) (8) 

In Eq.(^), we used the standard notation in which the lower indices refer to the space on 
which the matrices are acting and d 2 ± = Pd i2 P with P the flip operator. Eq.(|7|) can be 
proved by computing both sides and comparing them. 

Alternatively, eq.(^) can be written in a quadratic form. Let M 12 = J2ij-,ki ^ij-,kiEij ® E ki 
with, M ijlkl = -| (j^f-) SjkSu, and define 

r 12 = M 12 -M 21 -(MH-M*l) 
S12 = M 12 + M21- {M{\ + Ml\) 

Then, we have: 

{Vx f V 2 ) = r l2 V x V 2 + V 2 8vM. - V x s 12 V 2 - V 1 V 2 r 12 (9) 

Examples of this type of quadratic brackets first appeared in . Finally, we recall that the 
existence of a (i-matrix is equivalent to the property that the eigenvalues of the matrix V 
are Poisson commuting |13| . 



3 The N-soliton solution as a N-body problem. 

From the definition of the sine-Gordon field in terms of the r-functions, we see that only the 
eigenvalues Qi of the matrix V are important. We have 




In this section, we rewrite the soliton dynamics in terms of the variables Qi. This leads 
to a relativistic generalization of the Calogero model, first introduced by Ruijsenaars. We 
consider for definitness the evolution with respect to the light cone variable z + = x + 1. Let 
U be the matrix which diagonalizes V: 

V = U~ 1 QU (11) 
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where Q denotes the diagonal matrix diag(Qj). We define 

L = U/jU- 1 (12) 

where /i denotes the diagonal matrix diag( fii). As above, /Zj are the rapidities of the solitons. 
The matrix L plays the role of a Lax operator. Obviously, the quantities tr(L n ) = J2iLi $ 
are conserved during the evolution of the solitons. Moreover, they are in involution under 
the Poisson bracket (^). Finally, we have 

Proposition. The "time" evolution of L is given by a Lax equation 



L=[M,L], M = UU 



-i 



Here the dot means . Remarkably, L and M can be expressed in terms of the quantities 
Qi and Qi. 



L - = 2 rnnT> and M - = ^tH 1 -^ ( 13 ) 

Proof. We start from the relation 

fjV + Vfi = 2|e >< e| 

where |e > is the column vector with components = y/XiJTi. This relation remains true 
when we go to the basis where V is diagonal. Using the definitions eqs.([ll],|l2|), we get 

QL + LQ — 2|e >< e|, with \e>—U\e> 

Since Q is diagonal, we immediatly obtain from this relation 

„ ei e i 



U 



'J 



Next, we remark that = u.i\/~Xi) 

V = nV + Vu.= \e><e\ = U' 1 (Q + [Q, M]) U 

Multiplying this equation on the right by U^ 1 and by U on the left, and using the definition 
of |e), we get \e >< e\ = Q + [Q , M}. In components, this reads 

eiCj = Qi,5ij + (Qi - Qj)Mij 

If i — j, we find = yOi, and if i ^ j, we find the value of My in terms of Qi, Qi. This 
completely determines the matrix M since we can always ensure that M is antisymmetric 
by a suitable choice of basis. This ends the proof. 
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4 A canonical transformation. 



We now come to the description of the Poisson structure in terms of the variables Qi, Qi. 
Instead of the variables Qj, let us introduce the variables pi = Qi/Qi- The Lax matrix 
becomes 



hj ~ 2 Q^Q^ 

Comparing with eq.(D), we see that L has exactly the same form as V with the change of 
variables (pi, Gtj) — > (Q iy pi). This symmetry extends at the level of the symplectic structure. 
We have 

Proposition. The transformation (pi,cii) — > (Qi,Pi) is a symplectic transformation. 
Proof. We want to show that 

u = ^ A ^ + £ M^ \^ A ^ (14) 

i—l 0>i Pi j<-j \Pi Pj J Pi Pj 

s _Pi A S _Q± + y ( Igigj \ 5 _9i A 6 Qi 

i=l Pi Qi i<j \Qi Qj J Qi Qj 

Let us take as independent variables the pairs (Qi,pi), and expand both lines of eq.(|14"D in 
these variables. The equality of the two-forms then requires that: 

= (Hp__*p. + _p!i\ =Q (15 ) 
y Oj dpj aj dpi pf — p^ J 

<? (Qi_dpi_ _Q 1 dpj_ AQ.Qj \ 

\PidQi pjdQj Q\-Q}) 
Qj ddj _ pj dpj 
<k dQj pj dpi 

Let us first prove the third identity. By differentiating at pi fixed the matrix V written in 
either forms eq.(^[) or eq. ([□]), we obtain: 

dQ+ [Q, dU U' 1 ] =UdVlf- 1 = - (([/ a^da U~ l ) Q + Q{U a~ l da U~ v 

with a the diagonal matrix diag(aj). Using the fact that both Q and a are diagonal matrices, 
we derive that: 

dQi _ sr^ 77 ^ Q J rr-l 

O ~ ,J a- J ' 

Vi j=\ u j 

or equivalently, 



VVr tt- 1 ( ® k daj 



>ik 
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Introducing the matrices Aiij = UijU^ 1 and Ajk = ^~'§q~i we can rewrite this equation in 
matrix form: 

MA = Id (16) 

Similarly, using the expression of L in terms of pi and Qj, we find, by differentiating the 
relation L = UpU^ 1 at Qj fixed, that: 

TIM = Id (17) 

where IZ^j = ^|^- Comparing eqs. fllED anf ( |ITD proves the third relation in ([15]) since the 
inverse of A4. is unique. 

Let us now prove that Ay = 0. Expanding the first line of eq.(|14D in the variables (/ij, Qi) 
we have 

ij Qi /'./ i- j I 1 ' I 1 1 

The Poisson brackets {Qi,Qj} can be computed by inverting the symplectic form. We find 

{Qi,Qj} = (A^A'A- 1 



We know from the existence of the <i-matrix (|7|), that these Poisson brackets vanish. There- 
fore Ajj = 0. The relations Ay = are functional relations between the variables /ij, a$ at Qi 
fixed. Since the variables (/ij, cij) and (Qi, Pi) play a completely symmetric role, the conditions 
Aij = are also satisfied. This ends the proof that the transformation (/ij, aj) — > (Qi, pi) is 
symplectic. 

Explicitly, the Poisson brackets between (Qi, pi) are: 

{Qi , Qi} = o 
{pi , Qj} = QjPihj 

( AQjQj \ 

{p^pA = -{ofZ^Jp^ 

It will be convenient to introduce a new set of variables pi, which are conjugated to the 
coordinates Qi. We define them by: 

^ expw a(t^ (is) 

The Poisson brackets are then canonical: 



{Qi , Qj} = {Pi , Pj} = 

{Pi > Qj } Qj &ij 
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5 Hamiltonians and equations of motion. 

The result of the previous section has an immediate consequence. The Poisson brackets of 
the Lax matrix L, given in eq. (|T3|) , are quadratic with a dynamical R- matrix: 

{L x f L 2 } = R 12 L 1 L 2 + LaSisLx- L X S X2 L 2 - L X L 2 R 12 (19) 

where R\ 2 and S\ 2 are defined as r i2 and S12 in section 2, but with \Xi changed into Qi. This 
follows directly from the fact that the transformation (fii,a,i) to (Qi,pi) is symplectic, and 
the fact that the Lax matrix L is identical to the matrix V but with (/ij, a«) changed into 

(Qi,Pi) 

This implies that the function T(x) = det(l + xL) is a generating function of commuting 
Hamiltonians: 

{T(x) , T(y)} = 
Expanding it in power of x, we find: 

N 

T(x) = l + ^x p H p (20) 



P =i 

N 



— Qk 



p=l fci<---<fep 1 ki<kj \Qk + Qkj J 

As is well known, the existence of a i?-matrix ensures that all the flow generated by these 
Hamiltonians admit a Lax pair formulation. In particular, the Hamiltonian generating the 
evolution in the light cone coordinate z + is H + = tr(L): 

j=i j=i \m ^ k / 

The evolution in the other light cone coordinate z_ is generated by the inverse of the Lax 
matrix. Its Hamiltonian is i?_ = tr(L _1 ). Using tr(L _1 ) = H N _i/det(L), we find: 

Notice that one goes from H + to H_ by changing the sign of pj. 

The equations of motion follow from these Hamiltonians. For the light cone coordinate 
z + , we have 



Qi PiQi 



4QiQk 



where the dot still denotes ^f-. The flow is the z_ direction is very similar: one just 
formally changes tl 
by Ruijsenaars ||. 



formally changes the dot by the derivative These are the equations of motion written 
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To write explicitly the equations of motion for the (anti-) solitons and the breathers, we 
have to disentangle the reality conditions on (Qi,Pi). For the sine-Gordon field ( |10D and the 



Hamiltonians (|21"| , |22| ) to be real, the coordinates (Qi,Pi) have to come by pairs (j, j), with 
Qj = —Qj and Pj = Pj- The case j = j corresponds to a soliton or an antisoliton, and the 
case j 7^ j to a breather. Therefore, we introduce new coordinates qj by Qj = ie qj such that: 

Im q s = for s a soliton 
Im qj = it for s an antisoliton 
Qb — Qb f° r ^ a breather 

Similarly the momenta p s and pj are real and pb is complex with p^ — p b . In these coordinates, 
the Hamiltonians H± become: 



H ± = Y,e ±P > n^th(^) 



3 k^j 

The equations of motion for the flow generated by H + read: 

^sinh(gj-g fe ) 

By construction, these evolution equations are linearized in the variables (/Xj,aj). 
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